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Optical control of nanoscale objects has recently developed into a thriving field of research with 
far-reaching promises for precision measurements, fundamental quantum physics and studies on 
single-particle thermodynamics. Here, we demonstrate the optical manipulation of silicon nanorods 
in high vacuum. Initially, we sculpture these particles into a silicon substrate with a tailored 
geometry to facilitate their launch into high vacuum by laser-induced mechanical cleavage. We 
manipulate and trace their center-of-mass and rotational motion through the interaction with an 
intense intra-cavity field. Our experiments show optical forces on nanorotors three times stronger 
than on silicon nanospheres of the same mass. The optical torque experienced by the spinning 
rods will enable cooling of the rotational motion and torsional opto-mechanics in a dissipation-free 
environment. 


Nanoparticles often exhibit unique optical, mechani¬ 
cal, or electro-magnetic properties because of quantum 
effects in confined geometries and low dimensions [1]. 
Complementary to that, our present study is part of a 
long-term effort to control the quantum properties of the 
objects’ motion [2-4]. First experiments demonstrating 
de Broglie wave optics with macromolecules [5] were trig¬ 
gered by the question whether the superposition principle 
of quantum mechanics holds on all scales. They have led 
to the observation of quantum interference with masses 
beyond 10,000 amu [6]. An even higher mass regime, 
which might give insight to the quantum-classical transi¬ 
tion, can be reached with novel coherent manipulation 
schemes [7]. Models of a spontaneous localization of 
the wave function [8], and non-standard effects of grav¬ 
ity [9, 10], will become relevant for delocalized particles in 
the mass range of 10 5 — 10 8 amu, and above 10 10 amu, re¬ 
spectively. Tests of such models will require neutral, size- 
and shape-selected, cold, and slow nanoparticles that are 
mechanically isolated from their environment. This has 
motivated new experiments to launch and cool dielectric 
nanospheres in optical tweezers [11-13], ion traps [14], 
and in free-flight [15]. 

Here we extend this research to rod-shaped dielectrics 
of tailored geometry and anisotropic polarizability. Once 
their rotational motion can be controlled sufficiently well, 
these nanorods may be suitable for realizing torsional 
optomechanics [16-18]. Earlier experiments have shown 
that optomechanical torques can be exerted on nanorods 
in solution using the polarization or orbital angular mo¬ 
mentum of a light field [19-22]. The manipulation of 
nanoparticles in a dissipation-free environment, however, 
has remained challenging [23] . The coupling between ro¬ 
tational and motional degrees of freedom was recently 
demonstrated with optically trapped birefringent micro¬ 
spheres in a low pressure environment [24]. 
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FIG. 1. Launch and detection of size- and shape- 
controlled silicon nanorotors. An array of nanorods is 
carved into a silicon wafer. Focused laser pulses (1 — 3mJ 
pulse energy, 6 ns duration and 532 nm wavelength) hitting 
the rear side of the substrate cause stress-induced break¬ 
ing of the nanorods at predefined trenches on the front side 
and launch them into the high vacuum. Individual, freely 
moving, and rotating rods need to pass a 500 fim. diameter 
and 4 mm long aperture before entering a strongly pumped 
high-finesse infrared cavity (A = 1560 nm, 400 W intra-cavity 
power, 330,000 finesse, wq — 65 fim waist). Their velocities 
and rotation rates are deduced from the scattered light, which 
is collected by a multi-mode fiber placed perpendicular to the 
cavity axis and the direction of field polarisation. 


We have prepared periodic arrays consisting of more 
than a million silicon nanorods per mm 2 by dry-etching 
crystalline silicon wafers (see Appendix A). A scanning 
electron microscopy image of such an array is displayed 
in Figure 2 A. The individual nanorods exhibit a length of 
795 ±17 nm and a diameter of 108 ± 16 nm corresponding 
to a mass of (1.0 ±0.3) x 10 10 amu. The sample was posi¬ 
tioned in a chamber evacuated to 10 -6 Pa underneath an 
optical cavity (see Figure 1). The backside of the sam¬ 
ple was locally heated by a pulsed laser beam focused 
to 100 /im, which desorbs the rods by laser-induced ther¬ 
momechanical stress (LITHMOS)[15]. Small kerfs etched 
into the nanorod base define the desired break-off condi- 
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FIG. 2. SEM micrographs of the nanosculptured silicon rods. Array of nanorods with underetched kerfs (A) before 
and (B) after LITHMOS desorption. C & D) The kerfs define the break points where the rods crack. The actual rod length 
appears shortened under the oblique viewing angle. Nanorods with well-defined geometries can be launched, sent through the 
cavity and collected on a sample plate (E). Scale bars: (A & B) 5 /xm, (C) 1 fi m, (D) 400 nm and (E) 200 nm. 


tions (see Figure 2C & D; details in Appendix A). In Fig¬ 
ure 2A and B we show an electron micrograph of a sample 
spot before and after the LITHMOS pulses. It demon¬ 
strates that the rods can be reproducibly broken off at 
the tailored constrictions. Figure 2E depicts a close-up of 
the etched conical tip of a rod after launch and recapture. 

The standing light wave field of a high-finesse cavity 
allows us to track the translational and rotational mo¬ 
tion of the particles. It is optically pumped by a linearly 
polarized, distributed-feedback laser locked close to the 
cavity resonance. At the laser wavelength of 1560 nm 
silicon exhibits a high relative permittivity and minimal 
absorption. In a homogeneous field, the polarisability 
assumes a maximum value of aw/Aits^ = 6.4 x 10 9 A 3 
and a minimum value of a±/47r£o = 9.8 x 10 8 A 3 when 
the rods are oriented parallel and perpendicular to the 
field, respectively [25]. Even for rotating rods, the 
polarizability averaged over all possible rotation axes, 
Qo:|| + |aj_) /47 T£o = 2.8 x 10 9 A 3 , is larger than for 
a silicon nanosphere of the same mass, <a S ph/47T£o = 
1.4 x 10 9 A 3 . 

We can trace each nanorotor using the light it scat¬ 
ters into the direction perpendicular to both the cavity 


axis and the field polarisation. We collect this light in 
a 1 mm multimode fiber placed at a distance of 200 fi m 
from the cavity center. The detected intensity depends 
on the rod’s position in the standing wave and also on 
its orientation (see Appendix B). 

When a symmetric rotor enters the cavity and moves 
freely along the cavity axis, we expect a modulation of 
the scattering signal at two distinct frequencies: One is 
the translational frequency z/ trans = 2n^/A of the particle 
passing the standing-wave nodes with velocity v z , the 
second one is twice the rotation frequency, u rot = 2/ rot . 

Figure 3 A displays the normalized scatter¬ 
ing intensity of a freely rotating nanorod, 
Sn = (Is/Ic) /max (Ig/Ic), where Is is the mea¬ 
sured scattering signal and Ic the simultaneously 
recorded intra-cavity intensity. Panel C shows that it 
agrees well with the theoretical expectations for light 
scattering at dielectric needles, see Appendix B. We 
provide a full comparison of the measured signal in 
panel A and theory in Appendix C Figure SI. The corre¬ 
sponding Fourier spectrum, depicted in panel B, exhibits 
the distinct frequency contributions of translation and 
rotation. 
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FIG. 3. Tracing the nanorotor kinematics. The motion of every single nanorod can be deduced from the light it scatters 
while passing the cavity. A) The signal as a function of time conveys three distinct time scales corresponding to (i) the vertical 
transit through the cavity mode (Gaussian envelope), (ii) the transverse motion across the standing wave (low frequency 
z'trans), and (iii) the nanorod rotation (high frequency z/ ro t)- B) These frequency contributions can also be identified in the 
Fourier spectrum of the signal. The particle displayed here exhibits the vertical velocity v x — 11.5 d= 0.5 m/s, the on-axis 
velocity v z = 0.77 d= 0.05 m/s and the rotation rate / ro t = 2.15 ± 0.03 MHz. The geometric collimation of the incident particle 
trajectories permits an unambiguous distinction of the translational and rotational motion. C) The measured scattering signal 
(blue solid line) is well explained by a simple theoretical model (red dashed line), see Appendix B & Figure SI. 


For slow rods we could observe cavity-induced trans¬ 
lational channelling, i.e. one-dimensional trapping along 
an anti-node of the standing light wave. One such case is 
displayed in Figure 4, where the scattering signal (panel 
A) differs significantly from Figure 3. When averaged 
over the rotational period of the rod, the scattering signal 
(red dotted curve) does not drop to zero while the parti¬ 
cle is close to the center of the Gaussian beam. During 
this time the frequency related to the transverse motion 
of the rod is influenced significantly. In panel B we re¬ 
construct the particle trajectory through the cavity mode 
from the time evolution of the light scattering curve [15]. 
This is reproduced in a simulation of the rod’s motion 
under the influence of the cavity field (see Appendix C 
Figure S2). 

The optical channelling effect benefits from the geo¬ 
metrically enhanced induced dipole moment, due to the 
strong anisotropy of the rods [25]. For silicon nanorods 
rotating in the plane perpendicular to the cavity axis 


the orientation-averaged polarizability in the light field, 
(o!|| + aj_) /2, is enhanced by a factor of 2.7 in compari¬ 
son to silicon spheres of the same mass. In Figure 4 we 
observe an enhancement of the trapping potential by a 
factor of 2.6. 

In addition to the optical force on the particle’s center 
of mass, the cavity field exerts a torque on the rods, which 
tends to align them along the field polarisation axis [26] . 
Figure 5 displays the scattering curve of a slowly rotat¬ 
ing nanorod manipulated in both its center-of-mass mo¬ 
tion and its rotation. In analogy to Figure 4, we observe 
transverse channelling, but the rotation rate (black cir¬ 
cles) is now influenced significantly by the optical torque. 
This indicates that the motional and rotational degrees of 
freedom exchange energy via the optical potential. The 
influence of the optical potential on the rotation rate is 
most pronounced at around 12 /is. 

In conclusion, we have presented a method to tailor, 
launch, track, and manipulate high-mass silicon nanorods 
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FIG. 4. Transverse optical channelling: A) Scattering signal of a nanorod temporarily captured in an antinode of the 
standing-wave field. In contrast to Figure 3, the modulation envelope averaged over the fast rotation (red dotted line) no 
longer drops to zero between 5 — 14 /jl s. This indicates that the particle is channelled. B) Reconstruction of the center-of-mass 
trajectory of the rod. Channelling occurs at a trapping frequency around 470 kHz. A silicon sphere with the same volume 
would only be trapped at a frequency of 290 kHz [15]. 
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FIG. 5. Optically induced torque: Scattering signal (blue curve) of a channelled nanorod, whose average rotation rate 
(black circles, right scale) is deduced from the separation of two adjacent scattering maxima (red dots), i.e. half a rotation 
period. The cavity field couples the translational and the rotational motion. We observe that the particle first spins down to 
600kHz before it speeds up again to beyond 900kHz (see also Appendix C Figure S3). 


with well-defined geometry and high aspect ratio. Light 
scattering inside a high-finesse infrared cavity allows us 
to follow the translational and rotational motion of the 
nanorotors in real time. For some of the rods, cavity 
assisted ld-trapping and even rotational forces could be 
demonstrated. Our results are in good agreement with 
theoretical expectations and show that the rod-like shape 
enhances the interaction between the particles and the 
cavity field significantly, compared to silicon nanospheres 


of the same mass. Optomechanical trapping and cool¬ 
ing of the center-of-mass motion [4, 11-15] will benefit 
from aligning the rods along the axis of polarisation. 
Recent studies on single particle thermodynamics [27- 
29] may be extended to rotating systems. Our results 
represent a first step towards realizing torsional nano- 
optomechanics [16-18] and rotational cooling, which may 
become applicable to delicate biological nanomaterials 
such as the similarly shaped tobacco mosaic virus [30, 31] . 
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APPENDIX 

A. Nanorod sculpting 

The silicon nanorod arrays were fabricated from a single crystalline Si wafer by adapting a previously described dry 
etching method [32]. The 380 ym thick, <100> cut and p-doped wafers exhibit a resistivity of 1-10 Ocm. They were 
cleaned by sonication, first in acetone then in isopropyl alcohol (IPA), each for 5 min. They were then thoroughly 
rinsed with deionized water and dried in a stream of N 2 . The clean Si wafers were spin coated with MM A resist 
(Copolymer resist EL9, MicroChem) at 4000rpm for 60s, followed by baking at 180°C on a hotplate for 2min. 
PMMA resist (Polymer resist A2, MicroChem) was deposited on the MMA layer by spinning at 2000rpm for 60s, 
followed by baking at 180°C on a hotplate for 2 min. 

An array of 250 nm diameter circular dots with 1 y m spacing was written using a Raith 150 ultrahigh-resolution 
e-beam lithography system (Raith GmbH, Dortmund, Germany). The patterned wafer was developed in MIBK/IPA 
1:3 for 1 min, followed by rinsing with IPA for 20 s and drying in an N 2 stream. A nickel dot array was prepared 
by e-beam evaporating a 100 nm thick nickel layer at a base pressure of 10 _7 Torr with a rate of 1 A/s. Finally, the 
remaining resist was lifted-off in an Acetone/IPA 1:1 solution, washed with IPA and dried. The nickel nanodot arrays 
served as masks in the following dry etching. 

Vertical silicon nanowire arrays were fabricated by applying time-multiplexed reactive ion etching in an Inductively 
Coupled Plasma Deep Reactive Ion Etching machine (ICP-DRIE, PlasmaTherm SLR 770). In order to form nanopil¬ 
lars with a well defined breaking point, we have implemented a Bosch process, i.e. passivation followed by anisotropic 
etching. A variation of the ratio between the etching and the passivation times varies the scalloping and leads to 
different rod diameters. Based on that we set up a three-stage etching process: Six ’passivation-etching’ sequences 
with a time ratio of 1:1 allowed us to form 700 nm long nanorod segments. In order to create the breaking points, the 
ratio of the time windows was logarithmically changed to 0.7:1 during four further steps. Finally, a wider base was 
formed by changing the interval ratio linearly to 2:1 in six further steps. The silicon rods were cleaned by removing 
the nickel caps chemically. 


B. Scattering theory 

In order to compute the normalized scattered light intensity Sn , the rods are modelled as thin, homogeneous, 
dielectric needles of length L and diameter D. Adopting the scattering theory for dielectric needles [33] to a standing- 
wave situation, we find that the light intensity in the direction perpendicular to both the cavity axis e z and the field 
polarisation axis e x is proportional to 

oc k A D A L 2 f— — ^) \e y x u int | 2 [Si -f 2cos(2 kz)S+S- + 5 2 1 e -2 ^ +y 
J-C \£r T 1/ 

with k = 2tt/X the wave number, wq the cavity waist, and S± = sine [n • (e z db e y ) kL/2]. Here, we denote the center- 
of-mass position of the rod by (x, y , z), and the orientation of the rod is determined by the radial unit vector n. The 
internal field points in the direction ui nt = 2e x + (s r — l)(n • e x )n. 


C. Simulation 


Figures SI, S2 and S3 display the simulated dynamics of a L = 800 nm long and D = 100 nm thick silicon rod with 
mass M and moment of inertia 0 = ML 2 /12, which rotates in the plane perpendicular to the cavity z- axis. In this 
case the rod can be treated as a sub-wavelength particle at position z(t), and its orientation with respect to the field 
polarization x-axis is described by n = (cos </>(£), sin </>(£), 0). Given a constant intra-cavity field amplitude Fo, the 
rod’s motion is governed by the following classical equations of motion: 


z(t) = 


Elk 

4 M 


[aj_ + (ay — aj_) cos 2 cf\ sin {2kz) e 


-2 


(vxt) 


4>(t) 


_ 2 ( Vx 

— (a|| — aj_) sin (2 <fi) cos 2 (kz) e w o 


with cjo the cavity waist, v x the vertical velocity of the rod, and o^|| ? j_ the polarizability components, as given in the 
main text. Particle trajectories that do not pass through the cavity center, but slightly off-axis, can be accounted 
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for by decreasing the field amplitude Eq below the cavity value ^AIc/csq. The corresponding scattering intensity 
is obtained by evaluating the expression given in Appendix B along the simulated trajectory. Since the scattering 
signal is here proportional to Sn oc cos 2 (kz) exp (— 2v‘? c t 2 , the center-of-mass trajectory of a freely rotating rod 

can be reconstructed from the measured scattering signal by averaging over the fast rotation period. We ensure that 
we capture the full transit through the cavity mode by carrying out each simulation over the time interval between 
±10uq/v x . In order to compare this to the measured data, the scattering signal must be normalized, and the time 
offset of the simulation must be adjusted, accordingly. 



FIG. S 1. Normalized scattering signal of a freely rotating nanorod: Simulated signal as a function of time (red 
dashed line) in comparison to experimental data (blue solid line, see Figure 3). The following parameters and initial values 
were assumed: Eq — 4.15 x 10 6 V/m, v x = 11.5 m/s, 2 ( 0 ) = —89.85/fc, i(0) = 0.74 m/s, 0(0) = 0.1 rad, 0(O)/27t = 2.14 MHz. 



FIG. S 2. Simulation of a ld-channelled nanorod: (A) Normalized scattering signal of a ld-channelled particle. (B) 
Simulated center-of-mass trajectory along the cavity axis (blue solid line), compared to the trajectory reconstructed from the 
experimental data (red dots, see Figure 4). The following parameters and initial values were assumed: Eq — 8.2 x 10 6 V/m, 
v x — 11.3 m/s, 2 ( 0 ) — —91/fc, i(0) s== 0.28 m/s, 0(0) = —0.4rad, 0(O)/27r = 1.685 MHz. The scattering behaviour (C) and 
trajectory (D) of a sub-wavelength silicon nanosphere of the same mass is simulated for identical parameters. It illustrates that 
spherical particles are subjected to weaker optical forces compared to rods. 
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FIG. S 3. Simulation of translational and rotational manipulation of a rod: This simulation qualitatively resembles 
the case displayed in Figure 5 of the main text. The motional and rotational degrees of freedom couple via the optical potential, 
most dominantly after approximately 15 /xs. At this point, the channelled rod escapes the trapping potential of an anti-node 
and falls into the adjacent one, while the rotation rate temporarily slows down. The following parameters and initial values 
were assumed: E 0 = 8.0 X 10 6 V/m, v x = 7.94 m/s, z(0) — —99/k, i(0) — 0.28 m/s, 0(0) = —0.10098 rad, 0 (O)/ 27 t = 810 kHz. 
After transit, we find a 7% higher rotation rate, whereas the on-axis velocity is reduced by 34%. 
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